Instability of the negative mass Schwarzschild naked singularity 

Reinaldo J. Gleiser and Gustavo Dotti 

Facultad de Matemdtica, Astronomia y Fisica, Universidad Nacional de Cordoba, 
Ciudad Universitaria, (5000) Cordoba, Argentina* 

We study the negative mass Schwarzschild spacetime, which has a naked singular- 
ity, and show that it is perturbatively unstable. This is achieved by first introducing 
a modification of the well known Regge - Wheeler - Zerilli approach to black hole 
\Q • perturbations to allow for the presence of a "kinematic" singularity that arises for 

negative masses, and then exhibiting exact exponentially growing solutions to the lin- 
£vq . earized Einstein's equations. The perturbations are smooth everywhere and behave 

^! nicely around the singularity and at infinity. In particular, the first order variation 

of the scalar invariants can be made everywhere arbitrarily small as compared to the 
^ ' zeroth order terms. Our approach is also compared to a recent analysis that leads 

Q^ ■ to a different conclusion regarding the stability of the negative mass Schwarzschild 

spacetime. We also comment on the relevance of our results to the stability of more 
£\j \ general negative mass, nakedly singular spacetimes. 

>■ 
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I. INTRODUCTION: THE NEGATIVE MASS SCHWARZSCHILD SOLUTION 

The positive mass theorems of general relativity establish that, under plausible assump- 
tions, asymptotically flat solutions of the Einstein equations with physically acceptable 
matter sources cannot have negative total mass. Another well known feature of general 
relativity, compatible with the positivity theorems, is the possibility of the presence, or 
^ eventual formation, of singularities in the spacetime manifold. However, in the most rele- 

vant examples where this happens, the singularity is surrounded by a horizon, so that it is 
unobservable to an observer outside the horizon. That this should happen generically in all 
physically relevant situations is the essential content of the versions of the Cosmic Censorship 
hypothesis. On the other hand, from a strictly mathematical point of view, there is plenty of 
solutions of the Einstein equations that are asymptotically flat and present a negative total 
mass, but, in general, they are also nakedly singular, which would make them physically 
unacceptable. A simple and particularly interesting example of this type of solutions is the 
negative mass Schwarzschild spacetime. Nevertheless, for different reasons, these solutions 
are not expected to arise from physically acceptable regular initial data. One should notice, 
however, that what is physically acceptable may change as more data is available. What 
we have in mind here is the recent suggestion of the existence of "phantom matter" , that 
is some field that carries negative kinetic energy, as a possible explanation for the observed 
acceleration of the scale factor of the Universe [1]. 

A different, but equally important consideration regarding negative masses is the manner 
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in which they interact with positive masses. This gives rise to different types of instabilities 
that have been analyzed in the literature. We refer to [2] for a recent review. It is rather 
remarkable, however, that until recently [2], the stability of the negative mass Schwarzschild 
solution, that, in a loose sense, may be considered as representing the field of an isolated 
negative mass source, had not been considered in the literature, although the essential 
machinery for this analysis, at least at the linear level, had been available for some time. 
The analysis carried out in [2] leads to the conclusion that if certain boundary conditions, 
assumed suitable for the problem, are imposed, the negative mass Schwarzschild singularity 
should be considered as stable under linear perturbations. However, as discussed in Section 
VII, the marginally stable solution of Zerilli's equation used to prove stability in [2], is only 
a twice different iable functions of the radial coordinate. As a consequence, the associated 
perturbed metric is not continuous and thus it is not a solution of the linearized vacuum 
Einstein's equations (it actually satisfies the field equations of a singular distribution of 
sources). Given this situation, we undertake in this paper a new analysis of the linear 
stability of the negative mass Schwarzschild solution. This is based on a modification of 
the standard Regge- Wheeler- Zerilli formalism (reviewed in Section II) that simplifies and 
makes more transparent the analysis of the perturbative equations for negative mass. This 
is given in Section III. Next, we consider the existence of "zero modes", i.e., static solutions 
of the perturbation equations. We give the general solution, that for each mode depends on 
two integration constants, and prove that there is no choice of these constants such that the 
solution can be considered as a "small" perturbation. We give also a physical interpretation 
for these solutions that supports our conclusions. In Section VI, based on this modified form 
of the formalism, and this is the main result presented in this paper, we prove that under 
quite general boundary conditions, all scalar perturbative modes are unstable. Our results 
and possible consequences are finally summarized in Section VIII. 



II. PERTURBATION EQUATIONS FOR THE SPHERICALLY SYMMETRIC 

SCHWARZSCHILD SOLUTION. 

The vacuum solution of Einstein's equations for the exterior of a bounded spherically 
symmetric source, first obtained by Schwarzschild, may be written in the form, 

ds 2 = -(\-^M\ dt 2 +(\- 2 -M\ l dr 2 + r 2 {d9 2 + sm 2 (9)d<p 2 ) (1) 

where 8 and are standard coordinates on the unit sphere, and M is a real (positive or 
negative) constant, that can be identified with the Newtonian mass of the source. The 
range of r is limited by the singularities of the metric coefficients in (1). For M > the 
singularity at r = 2M corresponds to a regular horizon, and the solution may be extended 
to r < 1M . The resulting space time is known as a "black hole". The horizon hides a 
curvature singularity for r = 0, rendering the portion of space time outside the horizon 
globally hyperbolic. Similarly, for negative mass (M < 0) there is a singularity for r = 0, 
but in this case there is no horizon, the singularity is "naked" , and the resulting space time 
is not globally hyperbolic. 

The rather strange properties of the spacetime described by (1) led to question the pos- 
sibility of their physical realization. Since the solution (at least for r > 2M) is static, a 
first step in this direction, carried out by Regge and Wheeler [3], and completed later by 



Zerilli [4], was aimed at establishing their stability, through an analysis of the behavior of 
linear perturbations. The Einstein equations lead in this case to a set of coupled linear 
partial differential equations for the functions representing the perturbations. Because of 
the spherical symmetry of (1) this functions may be grouped into "even" (or scalar) and 
"odd" (or vector) sets that do not couple to each other, and are characterized by their 
dependence on the angular coordinates (0,(f)). Still, the resulting perturbative equations 
contain, in general, some redundant information that can be eliminated through a choice 
of "gauge". A particular choice, that keeps full generality, was selected in [3], [4] and is 
known as the " Regge- Wheeler gauge". Independently of the gauge issue, the linear stability 
of the Schwarzschild space time, well known for M > 0, was also established for vector 
perturbations for M < in [2]. For this reason we shall restrict our analysis to the scalar 
perturbations. In the Regge- Wheeler gauge, the scalar (even) perturbations, for each angu- 
lar mode (£,m), are fully described by four functions H (r,t), Hi(r,t), H 2 (r,t) and K(r,t), 
in terms of which the perturbed metric takes the form, 
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where e is an auxiliary parameter, useful in keeping track of the order of the perturbation. 
Yijn = Y^ m (9, 0) are standard spherical harmonics on the sphere, and, on account of the 
spherical symmetry, we may restrict tom = without loss of generality. The linearized Ein- 
stein equations for the metric (2) imply H (r,t) = H 2 (r,t), and a set of coupled differential 
equations for Hi, H 2 and K. This can be solved in terms of the Zerilli function ij)z(r, t), by 
the replacements, 
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provided ipz is a solution of the Zerilli equation 

d 2 i> z d 2 ip z 



where, 
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and r* is the "tortoise" coordinate, related to r by, 

£=(i-— V en 

dr \ r J 

In what follows, for simplicity, we choose the integration constant relating r and r* such 
that, 

r* = r + 2M [ln(r - 2M) - ln(2|M|)] (8) 

Therefore, for M < 0, we have r* = for r = 0. 

We notice that the relations (3) can be inverted, giving, 

, . . r(r-2M) / dK\ r T ^ 

te( '- " = (A + l)(Ar + 3M) {"> - r lF) + —l K < 9) 

Therefore, there is a one to one correspondence between %pz and the Regge- Wheeler gauge 
functions, Hq, Hi, H2 and K. 

For the stability analysis one considers a time dependence of the form, 

^ z {t,r) = e*"V(r) (10) 

which turns (5) into the Schrodinger-like form, 

d ^ z + V(r)ij = u 2 ij (11) 



dr 2 

Very loosely speaking, stability is attained if (11) has no solutions with negative u 2 . More 
rigorously, this statement must be supplemented with the imposition of appropriate bound- 
ary conditions on ip(r). The resulting boundary value problem for M > turns out to be 
quite different from that for M < 0. 

The case M > has been considered in detail in the literature (see, e.g. [5]). In this 
case we have —00 < r* < +00 and V(r) is positive definite in this range. Then, imposing 
boundary conditions that imply that the perturbed space time remains asymptotically flat, 
with a regular horizon, it can be shown that (ll)is turned into a self adjoint boundary value 
problem, that admits only solutions with u> 2 > 0, and, therefore, the spacetime is linearly 
stable. 

For M < the problem is quite different. The domain in this case is < r* < +00 
and, besides the singularity of V in (6)at r = r* = 0, there is another one at the interior 
point r = — 3M/A. Also, it is not at all clear what boundary conditions are physically 
adequate for r* = 0. Finally, Zerilli's auxiliary function ipz{t,r) is generically singular (at 
r = — 3M/A) for smooth metric perturbations. An attempt at an appropriate formulation of 
the problem was recently given in [2] . Unfortunately, for reasons explained below in Section 
VII, we consider that there are some difficulties in that analysis that require a revision of 
their conclusions. 



III. A MODIFIED REGGE- WHEELER-ZERILLI FORMALISM 

In view of the difficulties indicated in the previous Section, we present here a slight 
modification of the Regge- Wheeler-Zerilli formalism that, as we shall show, allows for a 



simple analysis of the stability problem for M < 0, with the intriguing result that these 
solutions are unstable under arbitrarily small perturbations. Following [2], we define \x = 
—2M, but consider, instead of ip, the function, 

C(t,r) = (2Ar-3^ z (t,r). (12) 

From (5) we easily find that ( satisfies the equation, 
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and shares with ipz{t,r) its uniqueness and gauge invariance properties. We may regard 
(14) as a map ^ := £ — > (if, iii, #2), a possible inverse map being 

The following results follow straightforwardly from the definition of \l/, eqn (14): 

Lemma 1: the metric perturbation given by ^(((t,r)) is C°° (for r e (0, 00)) if and only if 
((t,r) is C°° and satisfies the condition 
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If (13) is solved for initial data ((t = t ,r) and d(/dt(t = t ,r) both satisfying (16), the 
solution will satisfy (16) for all t where it is defined. In other words, (16) is preserved by 
the evolution equation (13) 

Proof: if the metric perturbation (K,Hi,H 2 ) is C°° then so is (, as follows from (15), and 
we can use Taylor's theorem in (14) to prove (16) is required to assure that (K, Hi, H 2 ) are 



smooth at r = 3p/(2A). The preservation of this constraint by the evolution equation (13) 
follows from the fact that the operator O defined in (13) sends functions satisfying (16) to 
functions satisfying (16) □ 

Since the coefficients in (13) depend only on r, we may separate variables by setting 
C(t,r) = exp(iu)t)x{r) , which gives 
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This ODE, obviously, does not have the Schrodinger like form of (11). Moreover, the 
coefficients of x an d dx/dr are singular at r = 3/i/(2A). However, the nature of this 
singularity is rather subtle, and, as shown below, all solutions of (17) are regular around 
r = 3/i/(2A), irrespective of the boundary conditions for r = and r = oo. 

Lemma 2: For any A,/i and oo, all solutions of (13) are regular at r = 3/i/(2A). 
Proof: Let p = r — 3/i/(2A), then expanding (17) around p = we get 
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If x is analytic around p = 0: 



then (18) reduces to 
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In our case ai = —2 so (20) splits into a system of three equations in the unknowns c , Ci 
and c 2 and an infinite system on the Ck,k > 2. Since the determinant of the finite system is 



&i a\ 

b 2 bi + a 2 2a 1 + 2 
&3 ^3 + h bi + 2a 2 
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(for all p, A and a;), and the infinite system can be solved iteratively in terms of c$, (20) has 
a solution for arbitrary values of c and C3, i.e., two linearly independent solutions (for more 
details on this special case of singular point see e.g. [9], pp. 205-215) □ 



IV. STATIC SOLUTIONS 

Static solutions correspond to uo = 0. In this case (17) is reduced to 
d 2 x (4Ar 2 + 6Ar/i - 3fi 2 ) d\ (4A 2 r 2 - 6\rfi + 3/^ 



dr 2 



(2Ar — 3/x) (r + fi)r dr (2Ar — 3/x) (r + /i) r 2 



x = o 



(22) 



A straightforward computation shows that this equation has, in general, two independent 
solutions, one regular and the other with logarithmic divergences for r = 0, and for r = — \x. 
The regular solution is a polynomial of order £ + 2, given by, 
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The general solution of (22) may be written in the form, 
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where C±, C 2 and C3 are constants. Notice that this expression contains effectively only 
two arbitrary constants, because a change in C3 amounts only to a change in C\. We also 
remark that, since all solution of (22) are regular for r > 0, we must have Q(r) > for 
r > 0, for otherwise a zero in Q(r) would introduce singularities in \ through the integral 
term in (25), which must be absent in the general solution of (22). 



A. Boundary conditions 

Equations (24,25) provide the general time independent solutions for the perturbative 
zero mode, but we still need to specify boundary conditions. The question is then if these 
solutions, perhaps for some particular choices of C\ and C 2 , are acceptable as perturbations 
of the negative mass Schwarzschild metric. As a possible criterion, we explore the change in 
a particular curvature invariant, which we choose to be the Kretschmann invariant /C, given 
by, 
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where R a bcd is the Riemann curvature tensor. Replacing (24)- (25) in the appropriate ex- 
pression for /C, linearized to first order in the perturbations, we find that, in the r ~ 
limit 
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We must set C 2 = to ensure that the first order term does not grow faster than the zero 
order term as r — > + . On the other hand, if we compute /C for large r, we obtain 
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then Ci must also be set to zero in order to keep the perturbation smaller than the 
zero order term at large r. We must, therefore, conclude that there is no (non trivial) 
choice of C\ and C 2 for which a static perturbation will remain small in the entire range of r. 



B. A physical interpretation for the zero modes 



Even though there are problems in interpreting the zero modes given by (23)- (25) as 
perturbations of the negative mass Schwarzschild metric, they are static solutions of the 
linearized equations and, therefore, their physical interpretation is of interest. One way 
in which the problem may be analyzed is to consider the Schwarzschild metric in Weyl 
coordinates. The general form of the static axially symmetric vacuum metrics in Weyl 
coordinates is, 



where U = U(p, z) 
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where a comma indicates partial derivative. It can be checked that (31) is the integrability 
condition for the system (32), so that k is obtained by quadratures if U, satisfying (31) is 
given. The latter is the Laplace equation in flat cylindrical coordinates and, as is well known, 
allows for the construction of static vacuum solutions of the Einstein equations by taking 
for U different solutions of the Laplace equation corresponding to localized sources. The 
Schwarzschild space time is recovered by taking for U the Newtonian potential corresponding 
to uniform line source of finite length, centered at the origin, and laying along the z-axis. 
The explicit forms for U and k in this case are, 
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where p± = \J p 2 + (z + M) 2 . The constant M is identified with the Schwarzschild mass 
and can be positive or negative. Suppose now that instead of Us we consider, e.g., 

U(p,z) = U s (p,z) + U 1 (p,z;m) (34) 

where Ui(p,z;m) is a solution of (31) that depends smoothly on a real parameter m such 
that, 

U 1 (p,z;0) = (35) 

This implies that U(p, z) given (34) satisfies also (31) and can be used to construct static 
vacuum solutions of the Einstein equations. But, in accordance with (35), these solutions, 
at least locally, can be made to approach the Schwarzschild metric arbitrarily closely, by 
taking m small enough. This may fail, and in general it does, when we approach the line 
source of Us- This line source corresponds to the regular horizon for M > and, therefore, a 
singularity in U\ destroys that regularity, and U\ cannot be considered globally as a "small" 
perturbation. Similarly, for M < 0, the line source corresponds to a naked singularity, and 
a singularity in U\ may modify its structure, and can make it even "harder". To illustrate 
this discussion we consider for U\ the form, 

Ui (p, z; m) = (36) 

' p z + z z 



This represents (in Weyl coordinates) a point source at the origin. The resulting equations 
for the modified k(p, z) may be solved explicitly, for any m, but we are only interested 
in the change to order m. Actually what we are interested in is the form of the metric 
in Schwarzschild, rather than Weyl coordinates. The coordinate change that accomplishes 
this is well known (see, e.g. [8], for details). Performing it we find that in Schwarzschild 
coordinates the modified metric, to first order in m, takes the form, 

2M 2m (r - 2M) 

g tt = -1 + 
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This metric may be expanded now in the Regge-Wheeler-Zerilli functions. For the metric in 
the form (37), in the notation of [3], the only non vanishing functions are H , H 2 , K, and 
G, corresponding to even (scalar) perturbations for even £, as the result of the symmetry 
of Ui(p,z;m) under z <-> —z. This is not in the Regge- Wheeler gauge, but, from gauge 
invariance [10], the function ip(r,t) is in this case given by, 

1>(r,t) = j^MX(r,t) 

, . 2r{r-2M)( dK (r - 3M) K\ „ / , . „.dG\ . . 
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Going back to (37), if we specialize to I = 2, we find 
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and the corresponding expression for x( r ) results, 
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which coincides with (23)-(25) for an appropriate choice of integration constants. Similar 
results hold for the other even values of £. Non vanishing contributions for odd £ may be 

obtained starting from a different choice of Ui(p,z;m) that makes it odd under z < ► —z, 

with analogous results regarding the corresponding zero modes. Therefore, the zero modes 
are the linearized expressions for the modifications that result in the Schwarzschild metric 
by the addition of different (static) multipoles. We stress again that these modifications ap- 
proach only locally the Schwarzschild metric, because they modify non trivially the structure 
of its singularity. We shall add more comments on this feature in Section VIII. 

V. TIME DEPENDENT SOLUTIONS 

Going back to (17), we consider now the general, time dependent case, with uj ^ 0. We 
are interested in solutions that are regular in < r < oo. Therefore, we need to consider 
the behavior of xi near the singular points r = 0, r = 3///(2A) and r = oo. First we find 
that near r = 0, any solution of (17), for any uj, is of the form, 

x(r) ~ C\r + C 2 r ln(r) + . . . (41) 

where C\ and C 2 are constants, and . . . indicates terms that vanish faster for r = 0. This 
implies that all solutions of (17) vanish for r = 0, but, in general, dx/dr diverges as ln(r). On 
account of our previous discussion regarding acceptability of the perturbations, this means 
that we need to impose the condition that this divergence is absent, and therefore, only the 
regular solution associated with C\ can be considered for the stability analysis. This fixes 
our boundary condition for r = 0. The behavior at r = oo is also simple. For large r, 
the general solution of (17) is of the form (55) [11]. We remind the reader that, in view of 
Lemma 2 in Section III, x( r ) is analytic at r = 3/x/(2A) 

From the assumed time dependence, positive values for uj 2 correspond to an oscillating, 
and therefore stable, behavior for the perturbations. The interesting range for instability is, 
therefore, u 2 < 0. On account of (13), the solutions for uj 2 < will grow unboundedly for 
large r, unless we can eliminate the exponentially divergent terms. In the following section 
we analyze this problem in detail. 
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VI. UNSTABLE SOLUTIONS 



Unstable modes, correspond to solutions of (17) with uj 2 = —k 2 , that satisfy the regularity 
condition for r = 0, and are exponentially decreasing for large r. Therefore, finding these 
solutions turns this equation into a boundary value problem, determining the appropriate 
value or values of k, if they exist. 

As a first step, we considered a "shooting" approach, and, using a high precision numerical 
method, we found, for £ = 2, a numerical solution for k ~ 4//i. Since this value was so close 
to an integer times l//x, we replaced k = 4///, £ = 2 in (17), and found that indeed it has a 
rather simple solution for this value of k. Repeating the procedure for £ = 3, we found again 
a solution for k = 20/ fi. This suggested the possibility that we could construct solutions for 
all £, by a simple generalization of the results for £ = 2, and £ = 3. This was effectively the 
case, and we found that, 



X(r) 



re 



-kr 



(r + fi) 



k/i 



C X + C 2 



{2Xy-3fi) 2 e 2ky 
c s (y + n) 



2kfi+l 



dy 



(42) 



-k 2 , provided 



(43) 



where Ci, C 2 , and C3 are arbitrary constants, is a solution of (17) with uj 2 
the constant k is given by 

, 2A(A + 1) (I - 1)1(1 + l)(l + 2) 

k = = 

3// 6/z 

and, therefore, kfi is an integer for all £. Notice that (42) effectively contains only two 
arbitrary constants, because any change in C3 can be compensated by a corresponding 
change in C\. 

It is easy to see that the integral in (42) contributes logarithmic divergences for r = 0, 
and an exponentially increasing factor for r — > 00. Therefore if we choose C 2 = 0, we find 
that 

((t,r) = dre- kr (r + fi) k>1 e kt (44) 

with k given by (43), is a solution of (17), representing the general expression for the regular, 
bounded, unstable modes. Moreover, this result implies that all scalar modes (i.e. for all 
X) of the Schwarzschild negative mass naked singularity are linearly unstable. 

As a check of the consistency of our results, we replace (44) in (12), and this in(3), and 
obtain 

(A +l)(r + */)*" 



K(t,r) = 
H^t) = -H 2 (t,r) 



3// 



■ exp[k(t — r)] 

A(A + 1)[2(A + l)r + 3fj] r(r + /i)^ 1 
V 



(45) 



exp[/c(t — r)] 



which shows explicitly the regularity of the solution of the perturbative equations in the 
range < r < +00. The perturbations of the curvature are also well behaved. Replacing in 
the Kretschmann scalar, we find that it is unchanged from the unperturbed value, to order 
e. Furthermore, a complete set of quadratic algebraic invariants of the Riemann tensor [6, 7] 
was computed to find that they are all invariant to first order in e. However, we have found 
non trivial first order contributions to the differential invariants 

180/i 2 (r + /i) U{£ -1)(£+1)(£ + 2) 



R, 



abcd:e 



R 



abcd;e 



Y f 



tm \ 



(46) 



x \r 



[r£(£ + 1) + 3/i] (r + ^ exp[k(t - r)] 
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and 

bcdief _ 90/i 2 (65/i 2 + 12(yr + 56r 2 ) 

^abcd;ef -ft ■ - ^ 

_ 30efe (r + ^exp^-r)] ^ (r)yaM) (4?) 

Q e (r) = (I 2 (/ - 1) (/ + 2) (Z + 1)V + 3/i/ (/ + 1) (/ 2 + / + 14) r 2 

+3// 2 (l7/ 2 + 17Z + 56)r + 180/i 3 ) (48) 

and we may expect non trivial contributions in some of the many other differential invariants. 
Note that, both for (46) and (47), the 0(e) term grows slower than the O(e ) term as r — »■ + 
limit, and also as r — ► oo. 



VII. COMPARISON WITH PREVIOUS RESULTS 

In [2], the Hilbert space 

H = L 2 {r*,dr*) (49) 

is introduced together with Zerilli's equation (11). Since the LHS of (11) is a self adjoint 
operator in 7i, it admits a basis of eigenf unctions that allow to evolve any initial data 
(ip , dtip)\t=o G 7i <S> 7i. It was also shown in [2] that (for £ = 2) there is a uj 1 = solution 
■0 O E 7t of (11). It is important to review at this point how ip is constructed. For negative 
M = — fi/2 the £ = 2 Schrodinger potential 

y = 2(l + g) ^(8r-W2) + |/.'(2r- < ./2) | (M) 

V r / r 3 ^r — ^y) 

is singular at r c = 3/i/4. The general solution of (11) for c<j = [2] 
r(3/x 3 — 6r 2 // + 4r 3 ) 



^o = C 3 - 



/i 3 (4r - 3/i) 

>(13/i 3 + 12r 2 /i - 24/i 2 r) r^/r 3 - 6r 2 /i + 4r 3 ) fr + fx 



+ C 4 



3/i 3 (4r - 3/i) /i 3 (4r - 3/i) 



log 



(51) 



lies outside Tt, however, for C3/C4 = log(7/3) — 4/9 = q , both ip and its first derivative 
vanish at r c [2], then ip and can be matched at r c to the trivial 00 = solution giving 

/ / n J fair) ,0 < r <r c 

W) = < (52) 

[0 ,r>r c 

Note that ip is an cu = solution in 7i, whose second derivative has a discontinuity at 
r = r c . The existence of such a solution signals a critical value of q = C3/C4 separating 
unstable (cj 2 < 0) from stable solutions /or £/ie Zerilli equation in 7i. 

In the previous section we exhibited a C°° perturbed metric which is an unstable solution of 
the linearized Einstein's vacuum equation around a negative mass Schwarzschild spacetime, 



13 

and showed that it satisfies appropriate boundary conditions, both at the singularity and 
at infinity. However, we have not given a procedure to evolve arbitrary perturbations. The 
reason why we have not adopted the Hilbert space (49) is that, although (49) is a natural 
choice in other contexts, generic linear perturbations (2) of the metric of a negative mass 
Schwarzschild spacetime lye outside H, since the integral 



{i>z^z)= / |V> 



)r A 



|2 



dr 



;i + /Vr) 



r(r + ^ } f H2 -r^\+-K 



'"' (53) 



;i + /Vr) 



6(r-3/i/4) V */ 3' 

diverges due to the pole at r = r c = 3///4 (here we have specialized to £ = 2 and used (9) to 
write the L 2 (r*,dr*) norm of the perturbation in terms of the perturbed metric elements, 
which may be arbitrary C 2 functions for r G (0, oo)). Even compactly supported metric 
perturbations are excluded as initial data in the formalism developed in [2], unless they 
meet the rather artificial requirement that, for every £ > 2, its £ — th harmonic projection 
satisfies 



' Q r 



Zjx 







2(£-l)(£+2) J 

The Zerilli function ipz is extremely useful because it reduces a complicated PDE system 
to a much simpler differential equation. However, care must be exercised when deciding 
what kind of solutions of this differential equation one is interested in, since ipz has, for 
the negative mass Schwarzschild spacetime, a built in singularity, as seen, e.g., from (9). 
As an example, it follows from (2) and (3), that the perturbed metric from (52) has a 
discontinuity at r = r c . In particular, it will not give a linearized solution of the vacuum 
Einstein's equations, in spite of being a solution of Zerilli's equation. This implies that, 
although (52) gives a marginally stable solution for the Zerilli problem in 7i = L 2 (r*, dr*), 
it does not prove the existence of a marginally stable solution of the linearized Einstein's 
equations we are dealing with, since only C 4 solutions to the Zerilli equation are relevant 
to the linearized gravity problem. Our choice of initial data set is that C|t=o an d d(/dt\ t =o 
be C°° functions subject to (16). As follows from Lemma 1 in Section III, this is equivalent 
to requiring that the metric perturbation be C°°. Equation (13) is a hyperbolic differential 
equation in (t,r) from which the evolution of given initial data for Hi, H 2 and K can be 
computed. The singularities in the coefficients of dx/dr and \i f° r 2Ar = 3/i were shown to 
be irrelevant in Lemma 2 of Section III. The hyperbolicity of (13), (or the regularity and 
smoothness of the light cones for the background metric), imply that smooth initial data of 
compact support given at, say, time t = to, should evolve into initial data also of compact 
support for t — ti, at least for sufficiently small t\ — to, then compactly supported smooth 
C\t=t and d(/dt\ t= t satisfying (16) evolves into a compactly supported C\t=ti 
The foregoing discussion shows that (45) can actually be excited without a fine tuning of 
the initial condition. We look at solutions of (13) corresponding to a time dependence of 
the form exp(iut), with uu arbitrary complex, without the logarithm singularity {C 2 = in 
(41)). A higher order expansion around r = 



X(r,u}) oc 



A(l + A) 3 2A(1 + A) 4 ^. 5 
V H y -— — -r^ + 0(r 5 



(54) 



3/x 2 9^ 3 

reveals that there is no term in r 2 , and only the coefficients of r 5 and higher order depend 
on A and on (even) powers of u. We have already shown in Section III (Lemma 2) that 
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all solutions of (13) for this type of time dependence are smooth at r = 3/x/(2A). We may 
similarly obtain an asymptotic expansion for these solutions for r — > oo. These can be 
written in the form, 



X(r,u) oc e 



[ik(r— /i ln(r+fi)] 



2iA(l + A) -3kfj. 2A(1 + A) + 3ik/j, 



2k\ 



4uj 2 r 



+ C(l/r 2 



(55) 



where k = ±uj, and the remaining coefficients depend only on even powers of u. For real uj, 
and large r, these solutions represent (modulated) incoming and outgoing waves, depending 
on the sign of k/u. Since (13) is linear, any linear superposition of solutions is also a solution. 
In particular, we may consider superpositions of the form, 



X(r,t) 



JC(uj)e lwt x(r : uj)duj 



(56) 



which provide everywhere regular solutions (including 2Ar = 3/x, for a wide variety of func- 
tions K,(u). An interesting choice is, 



K(u) = Ce- (a, -'" o)V/4 



(57) 



where C, ujq, and q are constants. Then, for large q, we may use a steepest descent estimation 
of (56) to find that for K(u) of the form (57), near r = we have, 



X (r,t) -K^^e™ * 



A(l + A) 3 2A(1 + A) 4 



3/i 2 



9/i 3 



r 4 + 0{r' 



(58) 



where K\ is a constant, and terms of order r 5 or higher are polynomial in t. This result 
implies that the amplitude near r = is strongly suppressed for large t. For large r, on the 
other hand, we find, 



X(r,t) :y. K 2 exp[iu u ± }exp[-(u ± ) 2 /q 2 ] 

2iA(l + A) - 3kfi 2A(1 + A) + Ziku 



r — 



2k\ 



+ 



4uj 2 7 



+ 0(l/r 2 



where K 2 is constant, and, 



M ±( r ! t) — t ± [r — /j, ln(r + /i)] 



(59) 



(60) 



Since this expansion is valid only for large r, we find that solutions corresponding to 
u + are strongly suppressed for all t, while those for U- represent an outgoing wave 
packet, roughly centered at t — [r — //ln(r + /i), moving along a null geodesic of the 
unperturbed spacetime. All these solutions represent the evolution of some smooth initial 
data that can be obtained from (56), for t — 0. Clearly, these superpositions display a 
"stable" behaviour, since they do not contain modes that grow exponentially in time. 
We notice, however, that we may combine any of these solutions with an arbitrarily 
small, but nonzero contribution from the unstable mode given by eqs. (44)- (45), to 
obtain an equally smooth initial data at t — 0, whose evolution displays the exponentially 
growing mode. We thus find a very large set of smooth initial data leading to an equally 
smooth "unstable" behaviour, showing that such data need not be particularly "fine tuned" . 
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Another issue discussed in [2] is that of the gravitational energy of the perturbation. Here 
again we find a number of subtleties, apparently peculiar to the negative mass case. The 
unstable solution (45) lies outside 7i in [2], and therefore outside the scope of ref [2]. To no 
surprise, its "Zerilli energy" 



E, 



dip z 
dt 


2 
+ 


dip z 



V\l/>2 



dr* 



(61) 



diverges. It is argued in [2] that Ez agrees with the notion of gravitational energy of a 
perturbation from second order perturbation theory, i.e., 



E, 



G 



ill 



— / G$rf?dL 



'(3), 



V = r 1/2 d/dt, ( = d/dt, 



(62) 



G ab being the second order correction to the Einstein tensor. However, from (45) we get, 
for the (£, m = 0) mode, 



G® 1 = exp(kt) 



A(rJ)(Y m ) 2 + B(r,£) 



dY, 



(£,0) 



00 



where A and B are C°° functions of r in (0, oo), falling off exponentially as r 
behaving as c/r 3 for r ~ 0. This implies that both 



(63) 



oo, and 



r 2 A 



-dr 



and 



r 2 B 



-dr 



'o / Jo J 

converge, and thus Eq is finite for the perturbation (45), showing that Ez ^ Eq in the 

negative mass case. 

We may try to find an alternative Hilbert space and energy concept associated to a Zerilli-like 

equation 

^ + 0(0) = O, = A(r)^ + B(r)^ + C(r). (64) 



dt 2 



Yet, is not hard to see that, for unstable spacetimes, all unstable modes will either be 
outside the Hilbert space, or have zero energy. To see this suppose we choose any hermitian 
product <, > in an appropriate function space under which O is self-adjoint: < 0i, (90 2 >=< 
(90i,02 >• Note that the hermitian product may well involve and r— derivatives of 0. 
Adding =< <f> t , (fi u + 0<t> > to its complex conjugated, and using the self adjointness of O, 
we arrive to a conserved "energy" 



dE d r ,, , ,_ , , 

-77 = J. [< <t>U <t>t > + < <f>, 0(f) >] 







(65) 



This is essentially a generalization of the argument used in [5] to obtain the conserved 
energy (61) above. We would like to use it now to prove that, if the underlying spacetime is 
unstable, the unstable modes will either be non normalizable (and then E in (65) undefined, 
as happens with the energy Ez for (45)), or will trivially vanish. In fact, suppose there is a 
solution 4>(t,r) = exp(nt)4> (r) , K > and assume O normalizable under <, >, so that the 
energy in (65) is well defined. Since 

E(t) = e 2Kt [k 2 < 0o, 0o > + < 0o, O0o >] 
does not depend on t, it has to vanish. 
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VIII. COMMENTS AND CONCLUSIONS 

In this article we have shown that even if we only allow for the most restrictive boundary 
conditions, both as r — > + and r — »■ oo, the perturbation equations for the negative mass 
Schwarzschild metric contain, for every angular mode, solutions that grow exponentially in 
time. These solutions represent everywhere regular perturbations that can be made initially 
globally arbitrarily small and, therefore, imply that the metric is unstable at the perturba- 
tive level, for all angular modes. An interesting technical point of our derivations was the 
introduction of a modification of the standard Regge-Wheeler-Zerilli method, that led to 
equation (17) and its rather unusual properties. Our analysis of the zero modes indicates 
that rather than perturbations they should be considered as linearized approximations for 
exact static metrics that locally, but not globally approach the negative mass Schwarzschild 
metric. Although asymptotically flat, they actually contain a naked singularity where the 
curvature diverges even faster than for the Schwarzschild metric. From the perspective of 
the present analysis, these spacetimes should also be considered unstable, since, if they are 
considered as part of the unperturbed metric, one would be allowed to relax the boundary 
conditions at r = 0, and as can be seen from the derivations in Section VI, this would lead to 
unstable solutions for any k = \J —uj 2 . Our results also suggest that, even though the vector 
(linearized) perturbations have been shown to be stable [2], one would expect that these 
would also become unstable through the coupling to the scalar modes. Since vector modes 
are related to rotation, this again suggests that rotation would not make the singularity 
more stable, although going through this in detail would require higher order perturbation 
theory, or, perhaps, some application of the Teukolsky equation. 

Finally, there is an obvious implication of our results to the cosmic censorship conjecture, 
since the linear instability would imply that the negative mass Schwarzschild cannot be the 
final stage of an evolving spacetime. 
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